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ABSTRACT 



We calculate the temperature profiles of (thin) accretion disks around rapidly 
rotating neutron stars (with low surface magnetic fields), taking into account 
the full effects of general relativity. We then consider a model for the spectrum 
of the X-ray emission from the disk, parameterized by the mass accretion rate, 
the color temperature and the rotation rate of the neutron star. We derive con- 
straints on these parameters for the X-ray source Cygnus X-2 using the estimates 
of the maximum temperature in the disk along with the disk and boundary layer 
luminosities, using the spectrum inferred from the EXOSAT data. Our calcula- 
tions suggest that the neutron star in Cygnus X-2 rotates close to the centrifugal 
mass-shed limit. Possible constraints on the neutron star equation of state are 
also discussed. 

Subject headings: X-rays:binaries-X-rays:spectra -stars:neutron -stars:rotation 
-Cygnus X-2 



1. Introduction 

The soft X-ray spectra of luminous low-mass X-ray binaries (LMXBs) are beheved 
to originate in geometrically thin accretion disks around neutron stars with weak surface 
magnetic fields (see for e.g. White 1995). An important parameter in modeling these 
spectra is the maximum value of the effective temperature in the accretion disk. The 
effective temperature profile in the disk can be estimated (assuming the disk to radiate from 
its surface like a blackbody) if one knows the accretion energy released in the disk. In a 
Newtonian treatment, the innermost region of an accretion disk surrounding a neutron star 
with weak magnetic field will extend rather close to the neutron star surface. The amount 
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of energy released in the disk will be one-half of the total accretion energy, the other half 
being released in the thin boundary layer between the disk's inner edge and the neutron 
star's surface. This then gives the disk effective temperature {Tes) varying with the radial 
distance (r) as Tes oc r~^/^ and the maximum effective temperature (T^^) will depend on 
the (nonrotating) neutron star mass (M) and radius (R) as T^"^ oc {MM/R^f/"^, where M 
is the steady state mass accretion rate. The value of (T^g^^) in the disk, in this approach, 
occurs at a radial distance 1.36 R. 

Mitsuda et al. (1984) parameterized the disk spectrum by the maximum temperature 
of the disk, using the above formalism and assuming the mass of the neutron star is equal 
to 1.4 Mq . These authors assumed that the inner parts of the disk do not contribute 
to the X-ray spectrum, and suggested a multi-color spectrum for the X-ray emission 
from the disk. It was shown by these authors, that the observed spectra of Sco X-1, 
1608-52, GX 349+2 and GX 5-1, obtained with the Tenma satellite, can be well fitted 
with the sum of a multi-color spectrum and a single blackbody spectrum (presumably 
coming from the boundary layer). White, Stella & Parmar (1988) (WSP) suggested that 
the simple blackbody accretion disk model should be modified to take into account the 
effects of electron scattering. Using EXOSAT observations, these authors compared the 
spectral properties of the persistent emission from a number of X-ray burst sources with 
various X-ray emission models. This work suggests that either the neutron star (in each 
system considered) rotates close to equilibrium with the Keplerian disk, or that most of the 
boundary layer emission is not represented by a blackbody spectrum. 

For accretion disks around compact objects, one possibility is that of the accretion 
disk not being Keplerian in nature. For e.g. Titarchuk, Lapidus & Muslimov (1998) have 
formulated a boundary problem in which the Keplerian accretion flow in the inner disk is 
smoothly adjusted to the neutron star rotation rate. The generality of such a formulation 
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permits application even to black holes, but only for certain assumed inner boundary 
conditions. These authors demonstrate that there exists a transition layer (having an extent 
of the order of the neutron star radius) in which the accretion flow is sub-Keplerian. An 
attractive feature of this formalism is that it allows super-Keplerian motion at the outer 
boundary of the transition layer, permitting the formation of a hot blob that ultimately 
bounces out to the magnetosphere. This formalism (Titarchuk & Osherovich 1999; 
Osherovich & Titarchuk 1999a; Osherovich & Titarchuk 1999b; Titarchuk, Osherovich & 
Kuznetsov 1999) therefore provides a mechanism for the production of high frequency 
quasi-periodic oscillations (QPOs) observed in the X-ray flux from several LMXBs. Such 
effects, when taken into account, can modify the Newtonian disk temperature profile 
(Chakrabarti & Titarchuk 1995). 

There are several other effects which will modify the Newtonian disk temperature 
profile, such as the effects of general relativity and of irradiation of the disk by the central 
neutron star. The wind mass loss from the disk and the residual magnetic field near the 
disk's inner edge may also play a part in modifying the effective temperature (Knigge 1999). 
Czerny, Czerny & Grindlay (1986) calculated LMXB disk spectra assuming that a disk 
radiates locally as a blackbody with the energy flux detemined by viscous forces, as well as 
irradiation by the boundary layer, and took into account relativistic effects, some of them 
in an approximate way. The possible effects of general relativity were also discussed by 
Hanawa (1989), using the Schwarzschild (nonrotating) metric, assuming that the neutron 
star radius is less than the radius of the innermost stable circular orbit (rin = 6GM/c^), 
which they identified as the disk inner boundary. The color temperature was assumed to 
be higher than the effective temperature by a factor of 1.5. It was found by Hanawa (1989) 
that the observations are consistent with a geometrically thin, optically thick accretion 
disk, whose inner edge is at r = rin, r being the Schwarzschild radial coordinate. 
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An important dynamical aspect of disk accretion on to a weakly magnetized neutron 
star is that the neutron star will get spun up to its equilibrium period, which is of the order 
of milliseconds (see Bhattacharya & van den Heuvel 1991, and refereces therein). The effect 
of rotation is to increase the equatorial radius of the neutron star, and also to relocate 
the innermost stable circular orbit (for a corotating disk) closer to the stellar surface (as 
compared to the Schwarzschild case). These effects will be substantial for rapid rotation 
rates in a fully general relativistic treatment that includes rotation. Therefore, for accreting 
neutron stars with low magnetic fields, the stellar radius can be greater or less than the 
radius of the innermost stable orbit, depending on the neutron star equation of state and 
the spacetime geometry. The effect of magnetic field will be to constrain the location of the 
inner-edge of the accretion disk to the magnetospheric (Alfven) radius. In such a case. Tin 
would lose the astrophysical relevance as discussed here. However, this will be so only if the 
magnetic field strength (B) is large. The problem addressed in this paper refer to LMXBs 
which contain old neutron stars which are believed to have undergone sufficient magnetic 
field decay (Bhattacharya & Datta 1996). Clearly, for low magnetic field case, a number of 
different disk geometries will be possible if general relativistic effects of rotation are taken 
into account. These structural differences influence the effective temperature profile and 
the conclusions derived by Czerny, Czerny & Grindlay (1986) and Hanawa (1989) are likely 
to be modified. 

In this paper, we attempt to highlight the effects brought in due to general relativity 
and rotation of the neutron star on the accretion disk temperature profile and then 
apply this to the particular case of the X-ray source Cygnus X-2. For simplicity 
(unlike Titarchuk, Lapidus & Muslimov 1998), we assume the accretion disk to be fully 
Keplerian, geometrically thin and optically thick. We first give a theoretical estimate of 
the modifications in T™"^ that would result if inclusion is made of the rotational effects of 
general relativity, and illustrate these by taking representative neutron star equations of 
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state. We then consider a model for the spectrum parameterized by the mass accretion rate, 
the color factor, and the rotation rate of the accreting neutron star (assumed to be weakly 
magnetized). We derive constraints on these parameters for the X-ray source Cygnus X-2, 
for which we take the estimates of T^^, the disk luminosity and boundary layer luminosity 
from the analysis of WSP. A conclusion of our work is that the neutron star in Cygnus X-2 
has a rapid spin rate close to the centrifugal mass-shed limit. 

The format of this paper is as follows. In Section 2, we discuss the rotational general 
relativistic effects on the disk temperature, using a formalism given by Page & Thorne 
(1974) and the disk irradiation by the neutron star. The theoretical predictions for the 
temperature profiles with these effects taken into account are presented in Section 3. 
Section 4 deals with comparison with observations, and its implications for parameters of 
our model for Cygnus X-2. A summary and discussions are presented in Section 5. 

2. The Effective Temperature of the Disk 
2.1. Effects of General Relativity and Rotation 

The effective temperature in the disk (assumed to be optically thick) is given by 



where a is the Stephan-Boltzmann constant and F is the X-ray energy flux per unit surface 
area. We use the formalism given by Page & Thorne (1974), who gave the following general 
relativistic expression for F emitted from the surface of an (geometrically thin and non- 
self-gravitating) accretion disk around a rotating black hole: 



(1) 



F(r) 
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where 



(3) 



Here Tin is the disk inner edge radius, E, I are the specific energy and specific angular 
momentum of a test particle in a Keplerian orbit and Qk is the Keplerian angular velocity 
at radial distance r. In our notation, a comma followed by a variable as subscript to a 
quantity, represents a derivative of the quantity with respect to the variable. Also, in this 
paper, we use the geometric units c — G — 1. 

For accreting neutron stars located within the disk inner edge, the situation is 
analogous to the black hole binary case, and the above formula, using a metric describing a 
rotating neutron star, can be applied directly for our purpose. However, unlike the black 
hole binary case, there can be situations for neutron star binaries where the neutron star 
radius exceeds the innermost stable circular orbit radius. In such situations, the boundary 
condition, assumed by Page & Thorne (1974), that the torque vanishes at the disk inner 
edge will not be strictly valid. Use of Eq. (1) will then be an approximation. This will 
affect the temperatures close to the disk inner edge, but not the T^^^ to any significant 
degree (see section 5 for discussion). 

In order to evaluate T^s using Eq. (1), we need to know the radial profiles of E, I and 
JIk- For this, we have to first compute the equilibrium sequences of neutron stars in rapid 
rotation. These can be calculated by noting that the space-time around a rotating neutron 
star can be described by the following metric (Cook, Shapiro & Teukolsky 1994): 



where the metric potentials 7, p, a, and the angular velocity (a;) of the stellar fiuid relative 
to the local inertial frame are all functions of the quasi-isotropic radial coordinate (f) and 



= gx^dx^dx^ (A, /3 = 0,1, 2, 3) 

= -e^+Pde + e2"(cif2 + f'^de^) + e^-Pf^sin'^e{d(j) - ujdtf , 



(4) 
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polar angle (9). r is related to the Schwarzschild radial coordinate (r) through the equation 

On the assumptions that the matter is a perfect fluid and that the space-time described 
by metric (4) is stationary, axisymmetric, asymptotically flat and reflection-symmetric 
(about the equatorial plane), the Einstein field equations reduce to three non-homogeneous, 
second-order, coupled differential equations (for 7, p and lS) and one ordinary differential 
equation (for a) in terms of e and P (the total energy density and the pressure of neutron 
star matter respectively) in addition to terms involving 7, p, uj and a (see Komatsu, Eriguchi 
& Hachisu 1989). We have solved these equations (self-consistently and numerically) to 
obtain 7, p, a;, a, P and Q (which is the angular velocity of the neutron star matter as 
measured by a distant observer) as functions of f and Q. The angular velocity enters in the 
equations through the rotation law (which must be specified) for the matter distribution. 
The equilibrium solutions so obtained can then be used to calculate bulk structure 
parameters such as gravitational mass M, equatorial radius i?, angular momentum J, etc. 
of the rotating neutron star. We assume that the neutron star rotates rigidly. Thus, f2 is 
constant for the stellar matter distribution, and is taken to be equal to O*, defined as the 
angular velocity of the neutron star as measured by a distant observer. 

Eq. (1) gives the effective disk temperature T^^ with respect to an observer comoving 
with the disk. Prom the observational viewpoint this temperature must be modified, taking 
into account the gravitational redshift and the rotational Doppler effect. In order to keep 
our analysis tractable, we use the expression given in Hanawa (1989) for this modification : 

1 + ^ = (1-— )-^/^ (5) 
r 

With this correction for (1 + 2;), we define a temperature relevant for observations (Tobs) as: 

^obs = 7— — 2^efr (6) 
1 + 2; 
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2.2. Computation of E, I and Qk 

For the work presented in this paper, we compute constant gravitational mass (M) 
equihbrium sequences for rigidly and rapidly rotating neutron stars using the formalism 
described above (see Datta, Thampan & Bombaci 1998 for details), keeping in mind the 
importance of the parameters M and for modeling the X-ray emission from LMXBs. 
These sequences are constructed starting from the static limit all the way upto the rotation 
rate corresponding to the centrifugal mass shed limit. The latter limit corresponds to 
the maximum fl^ {—^ms) for which centrifugal forces are able to balance the inward 
gravitational force. We now briefly describe how the quantities E, I and Qk are calculated; 
for details, the reader is referred to Thampan & Datta (1998). For a material particle in the 
gravitational field described by metric (4), we can write down (see for e.g. Misner, Thorne 
& Wheeler 1973) the equation of motion in the equatorial plane. These will be in terms 
of E, I, cu, f and the metric coefficients. The equation of motion in the radial direction 
defines the effective gravitational potential. The two conditions for orbits (circularity and 
extremum) at any r yield values for E and I as given by: 

I = (8) 

where v = {VL — uj)fe~P sinO is the physical velocity of the matter. The equation of motion in 
the azimuthal direction and that in the time direction yield the Keplerian angular velocity 
as 

Qk = e^''(^-)J/^+a;(r-) (9) 

[E — LUl) 
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2.3. Computation of Ebl and £'d 



We define the specific gravitational energy release due to the ingress of a material 
particle from infinity to the disk inner edge as Ej^, and that due to the particle spiralling in 
from the disk inner edge to the surface of the star as the boundary layer energy: E^l- For 
the case where the disk inner edge coincides with the stellar surface, Ebl is the difference 
in the energy of the particle in a Keplerian orbit at r = i? and that when it is at rest on 
the stellar surface. The exact expressions for E-ql and Ed are determined by the effective 
potential corresponding to any given space-time metric. For the Schwarzschild metric 
and the 'slow'-rotation Hartle-Thorne metric, the boundary layer to disk luminosity ratio 
has been calculated by Sunyaev & Shakura (1986) and Datta, Thampan & Wiita (1995) 
respectively. Calculations of Ebl and Ed corresponding to the metric (4) and used for the 
modehng in this paper are discussed in detail in Thampan & Datta (1998). 



For luminous LMXBs, there can be substantial irradiation of the disk surface by the 
radiation coming from the neutron star boundary layer. The radiation temperature at the 
surface of a disk irradiated by a central source is given by (King, Kolb & Burderi 1996) 



where rj is the efficiency of conversion of accreted rest mass to energy, f3 is the X-ray albedo, 
h is the half-thickness of the disk at r and n is given by the relation h (x r'^. For actual 
values of /3, h/r and n, needed for our computation here, we choose the same values (i.e., 
0.9, 0.2 and 9/7 respectively) as given in King, Kolb & Burderi (1996). Although the above 
equation is derived based on Newtonian considerations, corrections due to general relativity 
(including that of rapid rotation) will be manifested through the factor 77. We have made 



2.4. Disk Irradiation by the Neutron Star 




(10) 
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a general relativistic evaluation of rj for various neutron star rotating configurations, 
corresponding to realistic neutron star EOS models, as described in Thampan & Datta 
(1998). Since Tirr(r) oc r~^/^ and Tes{r) oc r~^/^, Tirr will dominate over T^s only at large 
distances. The net effective temperature of the disk will be given by (see Vrtilek et al. 1990) 

Tdisk(r) = (i;1r(r)+71^r(r))^/' (11) 

For the Cygnus X-2 modeling presented here, we find that Tirr does not play any significant 
role. However, since this quantity has consequences for the disk instability, we calculate it 
using Eq. (10), and illustrate it for the rotating neutron star models considered here. 



3. Results for the Disk Temperature Profile 

3.1. Neutron Star Equations of State 

The neutron star EOS is an important determining factor for the structure parameters 
of the star. A variety of neutron star EOS is available in the literature, ranging from very 
soft to very stiff models. For the purpose of our calculation, we have chosen four EOS 
models: (A) Pandharipande (1979) (hypcrons), (B) Baldo, Bombaci & Burgio (1997) (AV14 
+ 3bf), (C) Walecka (1974) and (D) Sahu, Basu & Datta (1993). Of these, model (A) is 
soft, (B) is intermediate in stiffness and (C) and (D) are stiff EOS. With this representative 
choice of EOS, the results of our calculations are expected to be of sufficient generality. 



3.2. The Results 

We have calculated the disk temperature profiles for rapidly rotating, constant 
gravitational mass sequences of neutron stars in general relativity. For our purpose here, 
we choose two values for the gravitational mass, namely, 1.4 Mq and 1.78 Mq , the former 
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being the canonical mass for neutron stars (as inferred from binary X-ray pulsar data), 
while the latter is the estimated mass for the neutron star in Cygnus X-2 (Orosz & Kuulkers 
1999). It may be noted with caution (Haberl & Titarchuk 1995), that this value is not 
confirmed from X-ray burst spectral analysis. We use the value of M = 1.78 Mq for the 
illustration of our results, and leave the issue for future confirmation. In order to make a 
comparison with observations of Cygnus X-2, we need to calculate the values of the £^bl 
and Eq, and as functions of the stellar rotation rate (O*), for the above chosen values 
of the gravitational mass (M). 

In Table 1, we list the values of the stellar rotation rate at centrifugal mass shed limit 
(Jims), the neutron star radius (R), the radius of the inner edge of the disk (nn), -Bbl, -Ed 
and the ratio Ebl/Ed, T^"^"" k T^f, and r'^y^^, for the two mentioned values of M 
and for the different EOS models. The last nine computed quantities are given for two 
values of neutron star rotation rate, namely, the static limit (Q* = 0) and the centrifugal 
mass shed limit (Q* = ^^ms)- -E'd and ErmBL are in specific units (i.e. units of rest energy 
moc^, of the accreted partilce). The temperatures are expressed in units of mI(^ x 10^ K 
(where Mn — M/10^^ g s~^). From this Table it may be seen that for a given neutron 
star gravitational mass (M): (a) Qras decreases for increasing stiffness of the EOS model, 
(b) R is greater for stiffer EOS. (c) The behaviour of rjn depends on whether rms > -R or 
Tins < R and hence appears non-monotonic. (d) E^l for the non-rotating configuration 
decreases with stiffness of the EOS. For a configuration rotating at the mass shed limit, 
Ebl is insignificant, (e) In the non-rotating limit, E^ remains roughly constant for varying 
stiffness of the EOS model. However, for the rapidly rotating case, the value of -Ed decreases 
with increasing stiffness, (f) The ratio Ebl/Ed in static limit is highest for the softest EOS 
model. For the rapidly rotating case, this ratio is uniformly insignificant, (g) and 
T^f^ decrease with increasing stiffness of the EOS models. However, these values exhibit 
non-monotic variation with (see Fig. 5 for the first parameter), (h) The rest of the 
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parameters, namely, and are non-monotonic with respect to the EOS stiffness 
parameter. 

In Fig. 1, we display the variation of R (the dashed curve) and (the continuous 
curve) with for M = 1.4 Mq for the four EOS models that we have chosen. Prom 
this figure it is seen that for a constant gravitational mass sequence, for both soft and 
intermediate EOS models, rin > R for slow rotation rates whereas, for rapid rotation rates 
Tin = R- In other words, for neutron stars spinning very rapidly, the inner edge of the disk 
will almost coincide with the stellar surface. It may be noted that for the stiff EOS models, 
this condition obtains even at slow rotation rates of the neutron star. 

It is instructive to make a comparison of the temperature profiles calculated using a 
Newtonian prescription with that obtained in a relativistic description using Schwarzschild 
metric. This is shown in Fig. 2, for the EOS model (B) and M = 1.4 M© (the trend is 
similar for all the EOS). The vertical axis in this figure is T^s (in this and all other figures, 
the temperatures are shown in units of mI(^) and the horizontal axis, the radial distance 
in km. This figure underlines the importance of general relativity in determining the 
accretion disk temperature profiles; the Schwarzschild result for is always less than the 
Newtonian result, and for the neutron star configuration considered here, the overestimate 
is almost 25%. For the sake of illustration, we also show the corresponding curve for a 
neutron star rotating at the mass shed limit (curve 4, Fig. 5). The disk inner edge is at the 
radius of the innermost stable circular orbit for all the cases. Note that the disk inner edge 
should be at R for Newtonian case; but we have taken rin = 6GM/ as assumed in Shapiro 
& Teukolsky (1983). 

The effect of neutron star rotation on the accretion disk temperature, treated general 
relativistically, is illustrated in Fig. 3a and 3b. Fig. 3a corresponds to the EOS model 
(B). The qualitative features of this graph are similar for the other EOS models, and are 
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not shown here. However, the temperature profiles exhibit a marked dependence on the 
EOS. This dependence is illustrated in Fig. 3b, which is done for a particular value of 
fl^ = flj^s. All these temperature profiles have been calculated for a neutron star mass equal 
to 1.4 Mq . The temperature profiles shown in Fig. 3a do not have a monotonic behavior 
with respect to Q*. This behavior is a composite of two underlying effects: (i) the energy 
flux emitted from the disk increases with Q^, and (ii) the nature of the dependence of Tin 
(where T^g vanishes : the boundary condition) on fl^ (see Fig. 1). This is more clearly 
brought out in Fig. 4, where we have plotted of versus for selected constant radial 
distances (indicated in six different panels) and EOS (B). At large radial distances, the 
value Teff is almost independent of the boundary condition; hence the temperature always 
increases with in Fig. 4 f. 

The variations of Ed, £^bl, the ratio £^bl/-E'D and T^^"^ with fl^ are displayed in Fig. 
5 for all EOS models considered here. All the plots correspond to M = 1.4 Mq . Unhke 
constant central density neutron star sequences (Thampan & Datta 1998), for the constant 
gravitational mass sequences, Ed does not have a general monotonic behaviour with 
jimax j^g^g ^ behaviour akin to that of Ed (because of the reasons mentioned earlier). £^bl 
decreases with Q*, slowly at first but rapidly as Jl* tends to Jims- The variation of Ebl/Ed 
with respect to fi* is similar to that of Ebl- 

We provide a comparison between the effective temperature (Eq. 1) and the irradiation 
temperature (Eq. 10), in Fig. 6. We have taken r] = E^l + Ed- Fig. 6a is for Q^. = 
while Fig. 6b is for a higher fl^ = 6420 rad s~^. The curves are for the gravitational mass 
corresponding to 1.4 Mq for the EOS model (B). The irradiation temperature becomes 
larger than the effective temperature at some large value of the radial distance, the ratio 
of the former to the latter becoming increasingly large beyond this distance. For Ebl 
small compared to Ed (as will be the case for a rapid neutron star spin rate), irradiation 
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effects in tfie inner disk region will not be significant. Defining the radial point where the 
irradiation temperature profile crosses the effective temperature profile as r = rcross and the 
corresponding temperature as Tcross) we display plots of rcross and Tcross with fl^ respectively 
in Figs. (7 a) and (7 b). It can be seen that rcross increases with Q*, just as Es does, and 
hence the irradiation effect decreases with increasing Q*. Therefore Tcross decreases with 
increasing Q^. 

In Fig. 8, we illustrate the disk temperature (Tdisk) profile for EOS model (B) 
corresponding to M = 1.4 Mq for various values of Q^. We illustrate the variation of Tdisk 
with fl^ at fixed radial points in the disk in Fig. 9. The effect of Tirr on T^isk can be noted 
in Fig. 9f. 

4. Compcirison with Observations: Implications for Cygnus X— 2 

The X-ray spectrum will have two contributions: one from the optically thick disk 
and the other from the boundary layer near the neutron star surface. The spectral shape 
of the disk emission depends on the accretion rate. For M « 10^^ g s~^, the opacity in 
the disk is dominated by free-free absorption and the spectrum will be a sum of blackbody 
spectra. The temperature of the local spectra (with respect to a co-moving observer) 
will be equal to the temperature Tef[{r) at that radius. The observer at large distance 
will see a temperature To]js{r), which includes the effect of gravitational redshift and 
Doppler broadening, as mentioned in Section 2. At higher accretion rates (M ^ 10^^ g s~^) 
the opacity will be dominated by Thomson scattering and the spectrum from the disk 
is that of a modified blackbody (Shakura & Sunyaev 1973). However, for still higher 
accretion rates Comptonization in the upper layer of the disk becomes important leading to 
saturation in the local spectrum to form a Wien peak. The emergent spectrum can then be 
described as a sum of blackbody emissions but at a different temperature than Tobs- The 
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spectral temperature seen by a distant observer is the color temperature Tcoi- In general 
Tco\ = fir)Tohs where the function / is called the color factor (or the spectral hardening 
factor), and it depends on the vertical structure of the disk. Shimura & Takahara (1988) 
calculated the color factor for various accretion rates and masses of the accreting compact 
object (black hole) and found that / (1.8-2.0) is nearly independent of accretion rate 
and radial distance, for M ~ M-^dd- These authors find that for accretion rate ~ 10% of 
the Eddington limit, / ~ 1.7. More recently, however, from analysis of the high-energy 
radiation from GRO J1655-40, a black-hole transient source observed by RXTE, Borozdin 
et al. (1999) obtain a value of / = 2.6, which is higher than previous estimates used in the 
literature. With this approximation for Tcoi, the spectrum from optically thick disks with 
high accretion rates can be represented as a sum of diluted blackbodies. The local flux at 
each radius is 

= 1^7rB,{fT,s) (12) 

where is the Planck function. For high accretion rates the boundary layer at the neutron 
star surface is expected to be optically thick and an additional single component blackbody 
spectrum should be observed. 

The £^X05'74T' observations of Cygnus X-2 (Hasinger et al. 1986) have been fitted to 
several models by WSP. One of the models is a blackbody emission upto the innermost 
stable circular orbit of the accretion disk and an additional blackbody spectrum to account 
for the boundary layer emission. The spectrum from such a disk is the sum of blackbody 
emission with a temperature profile 

T oc r-^/\l - {nJrY'^f^ (13) 

WSP have identified this temperature as the effective temperature which, as mentioned 
by them, is inconsistent since the accretion rate for Cygnus X-2 is high (M ^ 10^* g s~^). 
However, as mentioned above, identifying this temperature profile as the color temperature 
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makes the model consistent if the color factor is nearly independent of radius. Moreover, 
the inferred temperature profile (i.e. Tobs = T^oi/ f) is similar to the one developed in 
previous section. Therefore, in this paper we assume that the maximum of the best fit 
color temperature profile T^f^ is related to the maximum temperature T^f^ computed in 
previous section by (T™f^ P:i fT^f^). Shimura & Takahara (1988) suggested a value of 1.85 
for the factor /, for an assumed neutron star mass equal to 1.4 Mq and M = lOMEdd, where 
MEdd is the Eddington luminosity, with the mass to energy conversion efficiency taken as 
unity. 

For the source Cygnus X-2, the best spectral fit to the data is when T^f'^ — 1.8 x 10^ 
K, Ld = 2.1 X 10=^^ ergs s"^ and Lbl = 2.8 x 10=^^ ergs s"^ (WSP), where Ld is the disk 
luminosity and Lbl, the luminosity in the boundary layer. The distance to the source and 
the inclination angle to the source have been estimated by Orosz & Kuulkers (1998) to 
be 8 kpc and 60° respectively. Prom these values one can obtain, using the formalism 
described in section 2, the angular velocity of the neutron star (Q*) for a given neutron star 
mass, accretion rate (M), color factor / and equation of state. However, in order to make 
allowance for the uncertainties in the fitting procedure and in the value of z, and also those 
arising due to the simphcity of the model, we consider a range of acceptable values for T™f^, 
Ld and Lbl- In particular, we allow for deviations in T^f" and the best fit luminosities: 
we take two combinations of these, namely, (10%, 25%) and (20%, 50%), where the first 
number in parentheses corresponds to the error in T™f^ and the second to the error in the 
best fit luminosities. Note that we neglect the irradiation temperature here, as Tdisk ~ T'eff 
at the inner region of the disk (the region where the disk temperature reaches a maximum). 
The constraints on M, fl^, and / are calculated for two values of the mass of the neutron 
star in Cygnus X-2, namely, IAMq and 1.78M0. We obtain a range of consistent values 
for M, and / (and hence, allowed ranges of different quantities). The procedure is as 
follows. 



-18- 



As described in previous section, we can calculate the different quantities (£'d, -Bbl, 
T^f^, R, Tin, etc.) as a function of Q^,. Taking the observed (or fitted) values for T™f^, 
Lbl and (Lbl + -^'d) with the error bars, we have two limiting values for each of these 
quantities. We assume a particular value for each of / and M, from which we obtain the 
corresponding fitted values of T^f^, Ebl and {Ebl + -Ed) by the relations Ebl — Lbl/M, 
Ebl + Eb = {Lbl + Lb)/M and T^^"" = T^r/f- By interpolation, wc calculate two 
corresponding limiting Q^'s (i.e., the allowed range in fl^) for each fitted quantity. We take 
the common region of these three ranges, which is the net allowed range in Q*. We do this 
for Af 's in the range 10"^^ MqV-^ to 10~^ MQy-^ (which is reasonable for LMXB's) with 
logarithmic interval 0.0001, for a particular value of /. If for some M, there is no allowed 
Q*, then that value of M is not allowed. Thus we get the allowed range of M for a particular 
/. Next we repeat the whole procedure described above for various values of /, in the range 
1 to 10. If for some /, there is no allowed M, then that / is not allowed. Thus we get an 
allowed range of /. Taking the union of all the allowed ranges of M, we get the net allowed 
range of M (and similarly the net allowed range of Jl*) for a particular EOS, gravitational 
mass and a set of error bars. The allowed range of Ebl, Ed, R, etc. then easily follow, 
since their general variations with respect to are already known. The results of this 
exercise for various equations of state is shown in Table 2. Prom Table 2 we can read off the 
allowed range in /, — Q*/27r, i^in = QK,in/27r, R, r^^, M. For e.g. for gravitational mass 
M = 1.4 Mq , an assumed uncertainty of 25% in the luminosity, and 10% uncertainty in 
color temperature, these are respectively (1.37-2.39), (0.736-1.755) kHz, (0.745-1.787) kHz, 
(11.3-20.7) km, (16.0-28.3) km, (11.2-34.6) Medd- On relaxing the conditions on luminosity 
and color temperature to 50% and 20% respectively, the corresponding ranges change 
to (1.16-2.97), (0.719-1.743) kHz, (0.742-1.755) kHz, (10.7-20.7) km, (15.6-28.4) km, 
(5.8-42.4) MEdd- 

The EOS model (A) is the softest in the sample. The maximum mass of neutron 
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stars (at = Qms) corresponding to this EOS is 1.63 Mq. So the constraint results for 
Cygnus X-2, using this EOS are done only for M = 1.4 Mq . Since the luminosity in the 
boundary layer is about 10% of the disk luminosity, the neutron star is expected to be 
spinning close to the maximum allowed value. This is reflected in our results by the ratio 
of Q^/Qras ~ 0.95. In all these cases, the neutron star radius happens to be larger than the 
innermost stable circular orbit. Hence the radius of the inner edge of the disk coincides with 
the neutron star radius. Therefore, the angular velocity of the particles in the disk inner 
edge will be very nearly equal to that of the neutron star. This implies that the viscous 
torque in the disk inner edge will not be very significant, and the use of the Page & Thorne 
(1974) formalism will not introduce any gross error in the constraint estimates presented by 
us. 

5. Summciry and Discussion 

In this paper, we have calculated the temperature profiles of accretion disks around 
rapidly rotating and non-magnetized neutron stars, using a fully general relativistic 
formalism. The maximum temperature and its location in the disk are found to differ 
substantially from their values corresponding to the Schwarzschild space-time, depending on 
the rotation rate of the accreting neutron star. We have considered a model for the spectrum 
of the X-ray emission from the accretion disk, parameterized by the mass accretion rate, 
the color temperature, and the rotation rate of the star. We have compared the maximum 
effective temperature in the disk and the accretion luminosities (corresponding to the disk 
and the boundary layer) with the results of spectral fitting for the X-ray source Cygnus 
X-2 (WSP), and derived constraints on these parameters for the neutron star in this X-ray 
binary. The main conclusion of our analysis is that the neutron star in Cygnus X-2 has 
a rapid spin rate (close to the centrifugal mass shed value), and that the system has a 
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fairly large accretion rate (several times 10^^ g s~^). The low luminosity of the boundary 
layer compared to that of the disk for Cygnus X-2 is consistent with the above conclusion 
that the neutron star in this system has a rapid rotation rate. The low value of the ratio 
Lbl/Lj) justifies our assumption that the radiation pressure is negligible in the disk, so that 
the geometrically thin approximation for the disk is reasonable. According to Shimura & 
Takahara (1988), the spectrum from the disk can be represented as a multi-color blackbody 
only if M > O.lMedd- Our results for Cygnus X-2 are in accord with this. Interestingly, if 
we take the lower value 1.7 for the color index / (Shimura & Takahara 1988), we obtain 
a consistent set of results, except for the stiffest EOS model (D). This suggests that the 
comparatively lower values of / would disfavor stiff EOS for neutron star matter. However, 
if we take the value oi f — 2.6, as reported by Borozdin et al. (1999), one would require an 
EOS model that is stiffer than the stiffest used here, or a mass greater than M = 1.78 Mq (if 
one uses the narrower limits on the luminosity and color temperature). On the other hand, 
if one were to use the broader limits, the hardening factor / — 2.6 is disallowed only by the 
softest EOS model. 

We have assumed here that the magnetic field of the neutron star is weak, which 
implies that the radius of the last orbit of the accretion disk should be much greater than 
the Alfven radius (ra) (e.g., Shapiro & Teukolsky 1983), 



where M is the mass of the neutron star, //30 is the magnetic moment in units of 10^° G 



M = IAMq, the magnetic moment /^ao << 3.4 x 10~^ or the magnetic field in the surface 
should be less than 10^° G. So the conclusions presented by us will be vahd for the neutron 
star magnetic field upto a few times 10^ G. 

In our analysis, we have assumed that the the boundary layer between the disk and 




(14) 



cm' 



and Ta is in cm. The above condition implies that for R ^ 15 km., M/M^ 



edd 



20 and 
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the neutron star surface does not affect the inner regions of the disk. This will be a valid 
approximation when the boundary layer luminosity is smaller than the disk luminosity, and 
the boundary layer extent is small compared to the radius of the star. The flux received at 
earth from this region is 

FBL = (27ri?^)cos^(^) (15) 

U TT 

where Ai? is the width of the boundary layer, D — % kpc is the distance to the source, 
9 = 60° is the inclination angle and Tbl is the effective temperature. Spectral fitting gives 
a best fit value for Fbl ~ 4 x 10~^ ergs sce^^ cm~^ and Tbl = Tco1(bl)//bl = 2.88//bl keV, 
where /bl is the color factor for the boundary layer and T'coi(bl) is the color temperature of 
the boundary layer. Using these values, Ai? ^ 0.2 km, which is indeed smaller than 
R provided the boundary layer color factor /bl is close to unity. This is supported by the 
work of London, Taam & Howard (1986) and Ebisuzaki (1987), who obtain /bl ~ 1-5. 

A few comments regarding the validity of the Page & Thorne (1974) formalism for 
accreting neutron star binaries are in order here. Unlike for the case of black holes, neutron 
stars possess hard surface that could be located outside the marginally stable orbit. For 
neutron star binaries, this gives rise to a possiblity of the disk inner edge coinciding with 
the neutron star surface. We have assumed that the torque (and hence the flux of energy) 
vanishes at the disk inner edge even in cases where the latter touches the neutron star 
surface. In the case of rapid spin of the neutron star (as we infer for Cygnus X-2), the 
angular velocity of a particle in Keplerian orbit at disk inner edge will be close to the 
rotation rate of the neutron star. Therefore, the torque between the neutron star surface 
and the inner edge of the disk is expected to be neghgible. Independently of whether or 
not the neutron star spin is large. Page & Thorne (1974) argued that the error in the 
calculation of T^s will not be substantial outside a radial distance To, where Tq is given 
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by To — Tin = O.lrin- In our calculation, we find that (which is the most important 
region for the generation of X-rays) is greater than Tq by several kilometers for all the cases 
considered. 

The shortest time-scale of the system is given by the frequency in the innermost stable 
circular orbit (i^in, Table 2: column 5). A periodic oscillation in the system should be at a 
frequency lower than i/in (unless the model invoked to explain the temporal behavior predicts 
substantial power in the second harmonic, i.e., vqpo ~ 2z/in). The maximum frequency of 
the kHz quasi-periodic oscillation (QPO) observed for Cygnus X-2 is 1005 Hz (Wijnands 
et al. 1998). The stiffest EOS, model (D), will then be disfavored since i^in < 1 kHz for 
this model. Further, the neutron star mass estimate in Cygnus X-2 (f« 1.78Mq, Orosz & 
Kuulkers 1998) is not consistent with the soft EOS model (A). Our analysis, therefore, 
favors neutron star EOS model which are intermediate in the stiffness parameters. 

We have not attempted to model the observed temporal behavior of the source, and in 
particular, the QPO observations. Beat frequency models identify the peak separation of 
the two kHz QPO observed with the neutron star spin rate. For Cygnus X-2 the observed 
peak separation is Au — 346 ± 29 Hz (Wijnands et al. 1998) which is smaller than the 
typical rotation frequencies calculated here. However, a pure beat-frequency model has 
been called into question due to several observations. For instance, Ai/ has been observed 
to vary by about 40% for Sco X-1 (van der Khs et al. 1997) and the kHz QPO frequencies 
have been found to be correlated with the break frequency (^ 20 Hz) of the power 
spectrum density. An alternate model, where the QPOs are suggested to originate due to 
non-Keplerian motion of matter in the disk (Titarchuk & Osherovich 1999; Osherovich & 
Titarchuk 1999a; Osherovich & Titarchuk 1999b; Titarchuk, Osherovich & Kuznetsov 1999) 
have been proposed. These authors have also demonstrated the model by applying it to 
particular sources. Inclusion of this Newtonian model into the framework of the calculations 
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mentioned in this paper require a parallel formulation within the space-time geometry 
chosen herein. 

X-ray binaries hke Cygnus X-2 are beheved to be the progenitors of the miUisecond 
pulsars. Therefore, the discovery of a pulsar with a period fa 1 ms will strengthen the 
model presented in this paper, in terms of a rapidly rotating accreting neutron star. 
X-ray spectral analysis of Cygnus X-2 and similar sources using data from recent satellites 
(e.g. BeppoSAX, ASCA, Chandra) are required to provide further support to the model 
presented in this paper. 

We thank Paul J. Wiita for reading the manuscript and suggesting improvements in 
the presentation. 
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Table 1: Theoretically computed parameters: Centrifugal mass shed limit (f^ms), the neutron 
star radius {R), the disk inner edge radius (rin), specific gravitational energy release in the 
boundary layer (E'bl) and in the disk (Ed), their ratio Ebl/Ed, the maximum effective 
temperature {T^^"^), the radial location {r^"^) in the disk corresponding to T^'^^, T^f^ (see 
text) and the radial location (r^^) corresponding to this. These values are listed for two 
values of M for all EOS models considered here (except for EOS model (A), where the 
maximum neutron star mass is less than 1.78 M© , so only M — lA Mq is considered). The 



number following the letter E represents powers of 10. 



EOS 




(A) 




(B) 




(C) 




(D) 


M 




1.4 Mq 


1.4 Mq 


1.78 Mq 


1.4 Mq 


1.78 Mq 


1.4 Mq 


1.78 Mq 


(10^ rad s-i) 




11.026 


7.001 


8.219 


6.085 


6.808 


4.652 


5.088 


R 
(km) 


n = o 

r\ Q 


7.46 
11.44 


11.01 

1 C To 

15.72 


9.84 

1 c 1 n 

15.19 


12.28 
17.26 


12.32 

1 T OO 

17. 2o 


14.74 

on T A 

20.74 


15.76 

oi 1 

21.16 


rin 
(km) 


n = o 


12.40 
11.44 


12.41 
15.72 


15.81 
15.19 


12.41 
17.26 


15.75 
17.28 


14.74 
20.74 


15.79 
21.16 


Eel 


n = o 

17 17jjig 


0.275 
9.0E-5 


0.153 
5.0E-5 


0.262 
5.0E-5 


0.128 
4.0E-5 


0.185 
1.4E-4 


0.097 
1.4E-4 


0.136 
6.0E-5 


Eu 


n = o 

17 17jjig 


0.057 
0.073 


0.057 
0.057 


0.057 
0.071 


0.057 
0.053 


0.057 
0.064 


0.055 
0.045 


0.057 
0.054 


Ebl/Eb 


n = o 

17 — 17ms 


4.809 
l.OE-3 


2.673 
9.0E-4 


4.574 
7.0E-4 


2.248 
8.0E-4 


3.239 
2.0E-3 


1.779 

3.0E-3 


2.387 
l.OE-3 


(K (g s-i)-i/4) 


n = o 

17 — 17ms 


267.9 
320.2 


267.9 
261.7 


237.1 
277.1 


267.9 
246.3 


237.1 
255.6 


258.6 
215.5 


237.1 
221.7 


(km) 


n = Q 

17 17ms 


19.76 
16.14 


19.76 
21.64 


25.18 
21.42 


19.75 
23.68 


25.06 
24.05 


21.13 
28.39 


25.16 
29.21 


ymax/^1/4 

(K (g s-i)-i/4) 


n = o 

17 — 17ms 


224.8 
258.6 


224.8 
224.8 


197.1 
224.8 


224.8 
212.5 


197.1 
212.5 


218.6 
190.9 


197.1 
190.9 


^max 
obs 

(km) 


n = o 

17 — 17ms 


22.29 
18.70 


22.31 
23.69 


28.45 
24.58 


22.31 
25.60 


28.30 
26.90 


23.44 
30.14 


28.41 
31.72 
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Table 2: Observational constraints for various EOS : (A), (B), (C), (D). L and U stand for 
lower and upper limits. The parameters are f (color factor), z/* (frequency of the neutron 
star), i/in (frequency of the last orbit in the disk), R ( radius of the neutron star), r^'^ 
(radius where the effective temperature of the disk is maximum) and M (the accretion 
rate). The limits are for 25% uncertainty in luminosity and 10% uncertainty in the color 
temperature. Values in [ ] are for 50% uncertainty in luminosity and 20% uncertainty in 
the color temperature. For EOS (A), the mass of the neutron star cannot exceed 1.63 Mq 
hence the 1.7SMq solution is not presented. Medd is the Eddington accretion rate, which is 
1.4 X lO^^M/M© g s~^, where M is the neutron star mass. 



EOS 


M 




/ 






R 


„max 


M 




Me 






kHz 


kHz 


km 


km 


Medd 


(A) 


1.4 


L 


1.37[1.16] 


1.753[1.743] 


1.755(1.755] 


11.3(10.7] 


16.0(15.6] 


11.2(5.8] 






U 


1.99[2.56] 


1.755(1.755] 


1.787(1.944] 


11.4(11.4] 


16.1(16.1] 


22.9(27.5] 


(B) 


1.4 


L 


1.53[1.29] 


1.106(1.087] 


1.132(1.123] 


15.2(14.3] 


21.0(20.0] 


13.8(7.2] 






U 


2. 18 [2. 74] 


1.112(1.113] 


1.177(1.285] 


15.6(15.6] 


21.5(21.6] 


27.0(33.5] 


(C) 


1.4 


L 


1.57[1.33] 


0.964(0.945] 


0.975(0.971] 


16.8(15.6] 


23.1(21.7] 


14.9(7.7] 






U 


2.24[2.81] 


0.968(0.968] 


1.015(1.134] 


17.2(17.2] 


23.6(23.7] 


29.3(36.5] 


(D) 


1.4 


L 


1.67[1.42] 


0.736(0.719] 


0.745(0.742] 


20.1(18.6] 


27.6(25.7] 


17.5(9.1] 






U 


2.38[2.97] 


0.740(0.740] 


0.779(0.876] 


20.7(20.7] 


28.3(28.4] 


34.6(42.4] 


(B) 


1.78 


L 


1.58(1.33] 


1.303(1.292] 


1.322(1.315] 


14.8(14.2] 


21.2(20.7] 


8.9(4.7] 






U 


2.28[2.91] 


1.307(1.307] 


1.361(1.462] 


15.1(15.1] 


21.4(21.4] 


17.2(21.4] 


(C) 


1.78 


L 


1.65[1.39] 


1.081(1.067] 


1.086(1.085] 


17.1(16.2] 


23.8(23.0] 


9.8(5.1] 






U 


2.39[3.01] 


1.083(1.083] 


1.109(1.209] 


17.3(17.3] 


24.0(24.1] 


19.3(24.0] 


(D) 


1.78 


L 


1.74 [1.47] 


0.806(0.791] 


0.817(0.813] 


20.6(19.2] 


28.6(27.1] 


11.4(6.0] 






i: 


2. 50 [3. 15] 


0.809 [0.809] 


0.818(0.938] 


21.1(21.1] 


29.1(29.2] 


22.2(27.7] 
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Fig. 1. — Disk inner edge radius (rin, curve 1) and neutron star radius (i?, curve 2), as 
functions of neutron star angular velocity (Q*) for various EOS models. The curves are for 
a fixed gravitational mass (M = 1.4 M© ) of the neutron star. 
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Fig. 2. — General relativistic corrections to Newtonian temperature profiles for EOS model 
(B) and neutron star gravitational mass M — lA Mq . Curve (1) corresponds to the 
Newtonian case, curve (2) to the Schwarzschild case and curve (3) to a neutron star rotating 
at the centrifugal mass shed limit, calculated using the metric (4). For curve (1) it is assumed 
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Fig. 3. — Temperature profiles incorporating the effects of rotation of tlie neutron star. Tlie 
plots correspond to (a) EOS model (B) and an assumed neutron star mass of M = 1.4 M© for 
rotation rates: = (curve 1), il* = 3.647 x 10^ rad s'^ (curve 2), = 6.420 x 10^ rad s~^ 
(curve 3), = 7.001 x 10^ rad = ilms (curve 4) (b) the same assumed mass and fl^ = ilms 
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Fig. 4. — Plot of Tes versus for chosen constant radial distances for fixed neutron star 
mass M = 1.4 M© and EOS (B). The plots correspond to: (a) r = 13 km, (b) r = 18 km, 
(c) r = 35 km, (d) r = 100 km, (e) r = 2000 km, (f ) r = 5000 km. 
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Fig. 5. — The variations of the £^d, Ebl, Ebl/Ed and with for a chosen neutron 
star mass value of 1.4 M© for the four EOS models. The curves have the same significance 
as Fig. 3b. 




Fig. 6. — Comparison between the radial profiles of Tejf (curve 1) and Ti„ (curve 2), calculated 
for 77 = E-Bh + Ed, (3 — 0.9, h/r — 0.2 and n — 9/7 in Eq. (10), for two values of neutron 
star spin rates: (a) = and (b) Jl* = 6.420 x 10^ rad s~^. The curves are for a neutron 
star configuration having M = 1.4 Mq , described by EOS model (B). The temperatures are 
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Fig. 7. — Plots for: (a) rcross versus Jl* and (b) Tcross versus O*. These are for fixed neutron 
star gravitational mass of M = 1.4 Mq and for the different EOS models as in Fig. 3b. Here 
Tirr is calculated for rj = Ebl + ^d, P = 0.9, h/r = 0.2 and n = 9/7. 
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Fig. 8. — The disk temperature (Tdisk) profiles for a M = 1.4 Mq neutron star corresponding 
to EOS model (B) having various rotation rates as in Fig. 3a. These curves are obtained for 
77 = Ebl, and the same values of P, h/r and n as in Fig. 6. 
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Fig. 9. — Plots of Tdisk versus at various chosen radial distances: (a) r — 13 km, (b) 
r = 35 km, (c) r = 100 km, (d) r = 200 km, (e) r = 250 km, (f ) r = 1000 km. These are for 
EOS model (B), an assumed gravitational mass value of 1.4 Mq , and the same values of 77, 
/3, h/r and n as in Fig. 7. 



